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ORIGINAL ARTICLE 

Multiscale Design of Cell-Type-Specific Pharmacokinetic/ 
Pharmacodynamic Models for Personalized Medicine: 
Application toTemozolomide in Brain Tumors 

A Ballesta 1 , Q Zhou 2 , X Zhang 1 , H Lv 1 and JM Gallo 1 

Optimizing anticancer therapeutics needs to account for variable drug responses in heterogeneous cell populations within 
the tumor as well as in organs of toxicity. To address cell heterogeneity, we propose a multiscale modeling approach— from 
in vitro to preclinical and clinical studies— to develop cell-type-specific pharmacokinetic-pharmacodynamic (PK-PD) models. 
A physiologically based mechanistic modeling approach integrating data from aqueous solutions, U87 glioma cells, mice, and 
cancer patients was utilized to characterize the brain disposition of temozolomide (TMZ), the cornerstone of chemotherapy 
against glioblastoma multiforme.The final model represented intracellular normal brain and brain tumor compartments in which 
TMZ pH-dependent conversion to the DNA-alkylating species leads to the formation of DNA adducts that serve as an entry point 
for a PD model. This multiscale protocol can be extended to account for TMZ PK-PD in different cell populations, thus providing 
a critical tool to personalize TMZ-based chemotherapy on a cell-type-specific basis. 
CPT Pharmacometrics Syst Pharmacol. (2014) 3, e1 12; doi:1 0.1 038/psp.2014.9; published online 30 April 2014 



Personalized cancer therapy is a new horizon that offers to 
dramatically improve patient responses and early attempts 
to tailor chemotherapy are either rooted in a patient's 
genomic signature 1-3 or based on systems pharmacologi- 
cal approaches involving protein networks and pharmaco- 
kinetic-pharmacodynamic (PK-PD) models. 4-7 Regardless 
of whether genomic- or protein-centric methods are under- 
taken, an account of cell heterogeneity is in order. Indeed, 
recent analyses have revealed heterogeneity both within the 
tumor region — a collection of cell types, such as cancer cells, 
tumor stem cells, endothelial cells, and immune cells 8 — and 
within organs of toxicity that may include different healthy cell 
types; for instance, bone marrow contains hematopoietic pro- 
genitors of different maturities. 9 Each cell population may dis- 
play specific PK-PD characteristics and drug sensitivity that 
will impact the overall response. 

To address cell heterogeneity, a systems pharmacological 
approach utilizing network PK-PD models is proposed as a 
comprehensive mean to design patient chemotherapy. We 
considered the importance of constructing physiologically 
based models with intracellular compartments as this pro- 
vides relevant PK-PD features that could differentiate drug 
action in each cell type. However, the experimental assess- 
ment of such cell-type-specific intracellular PK-PD models 
remains challenging in a preclinical setting and to a greater 
extent in cancer patients. In particular, although one can 
readily appreciate that intracellular drug concentrations are 
the final input to drug action or PDs, there has been no tan- 
gible means to obtain this information in a whole animal. In 
this context, we developed a multiscale — in vitro to in vivo— 
modeling approach to bridge this gap. 

The multiscale modeling protocol was developed for temo- 
zolomide (TMZ), the cornerstone of chemotherapy against 



glioblastoma multiforme (GBM), the most frequent primary 
brain tumor in adults. 1011 GBM therapeutic management 
involving surgery, radiation, and chemotherapy has shown 
suboptimal advances as patient median survival remains at 
12 to 14 months and may benefit from personalized TMZ- 
based chemotherapy. TMZ is a prodrug that spontaneously 
converts into its metabolite 5-(3-methyltriazen-1-yl)imidaz- 
ole-4-carboxamide (MTIC), which is subsequently degraded 
into 4-amino-5-imidazole-carboxamide (AIC) — an inactive 
metabolite — and a methyldiazonium cation, the DNA-meth- 
ylating species. 12-14 These reactions are pH-dependent as 
TMZ and MTIC degradation rates increase and decrease, 
respectively, with pH values. 15 The methyldiazonium cation 
creates DNA adducts — a marker of TMZ PD — that trigger 
DNA damage responses involving the p53 pathway, DNA 
repair mechanisms, and potentially cell-cycle arrest and cell 
apoptosis. 14 ' 16 To initiate the development of cell-type-specific 
PK-PD models, we considered a homogeneous population of 
glioma cells for the brain tumor that was cast into intracellular 
PK models to ultimately specify the brain disposition of TMZ 
and its metabolites. 

RESULTS 

The multiscale approach to design a model of TMZ intracel- 
lular brain disposition progressed through the following steps: 
(i) TMZ pH-dependent conversion study in solutions, (ii) 
in vitro JUZ conversion and transport studies in U87 glioma 
cells, (iii) TMZ PK studies in mice, and (iv) TMZ PK studies in 
brain tumor patients. They are presented in sequence below. 

TMZ pH-dependent conversion 

The model of TMZ pH-dependent conversion was derived from 
the literature 12-14 and represents a stoichiometric degradation 
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of TMZ to its metabolite MTIC — rate constant k T — that sub- 
sequently fragments into the methyldiazonium cation and 
the inactive species AIC via rate constant k M (Figure 1a, 
equations in Supplementary Data). The rates k T and k M 
are highly pH-dependent since TMZ is stable at acidic pH but 
decomposes to MTIC at more alkaline pH values, whereas 
MTIC degradation rate is lower at large pH values and 
increases as pH decreases. 1215 They were modeled as: 



/c T (pH) = /c T0 -e A TPH 

MpH) = /c M o-e-^ H 

The instability of the methyldiazonium cation, which is a 
highly reactive species was represented by a degradation of 
rate constant k cat set to 6,000 hr 1 that is equivalent to a half- 
life of 0.4 s. 17 Parameters k J0 , X v k m , and X M were estimated 
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Figure 1 A multiscale approach for mechanistic modeling of TMZ PK. (a) TMZ pH-dependent conversion in solution. Each molecule of 
TMZ produces one molecule of MTIC that further degrades into one molecule of AIC, an inactive metabolite, and one molecule of the 
methyldiazonium cation, which is the methylating agent, (b) TMZ PK and simplified PD model in U87 glioma cells. See Results section for the 
definition of mathematical variables and parameters, (c) Model of TMZ brain disposition in tumor-bearing mice and brain tumor patients. The 
blood compartment is a forcing function, whereas the normal brain and brain tumor compartments are represented by physiologically based 
models. See Results section for the definition of mathematical variables and parameters, (d) Prospective cell-type-specific model of TMZ 
PK-PD.The brain disposition model can be extended to account for cell heterogeneity both in the tumor region and in the bone marrow, the 
main toxicity target of TMZ. The blood and interstitial fluid compartments are directly inferred from the brain disposition model, whereas each 
cell type is represented by a model of TMZ intracellular PK-PD fitted to in vitro data performed in the corresponding cell population. MTIC, 
metabolite 5-(3-methyltriazen-1-yl)imidazole-4-carboxamide; PD, pharmacodynamics; PK, pharmacokinetics; TMZ, temozolomide. 
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utilizing available TMZ, MTIC, and AlC concentration-time 
data in buffer solutions at four pH values (1, 6.8, 9.1, and 
13) with an initial TMZ concentration of 1 mmol/l 15 (Table 1). 
The best-fit model closely agreed with the observed data 
(Figure 2a-d). 

For TMZ exposures ranging from 2 to 96 h, the percent- 
age of activated TMZ-methylating cation had the form of an 
increasing steep S-shape function as pH increased, suggest- 
ing that a small pH variation may lead to a significant change 
in the TMZ activation rate (Figure 2e). Further, the model 
predicted that TMZ activation rate was close to zero at pH 
= 4 for any exposure duration and increased with pH values 
(Figure 2a-d). Interestingly, although the rate of MTIC frag- 
mentation decreases with pH, the reaction does not prevent 
TMZ activation at high pH. 

In wfroTMZ PK 

The second step of the multiscale approach analyzed TMZ 
metabolism, membrane transport and simplified PD in U87 
glioma cell culture. Two physiological compartments — extra- 
cellular and intracellular — were considered in which TMZ 
pH-dependent activation and MTIC subsequent degradation 



were represented by the above-described model (Figure 1 b). 
Because TMZ is highly lipophilic and constitutes a poor sub- 
strate of ABC transporters, 18 transport of the parent drug and 
AlC between the extracellular and intracellular compartments 
were modeled as passive diffusion. As MTIC displays limited 
ability to cross cell membranes and as the methyldiazonium 
cation is a highly reactive species, their transport between 
the two compartments were not considered. 19 Concerning 
TMZ PD, the methyldiazonium cation is the sole species able 
to form DNA adducts. 12 

Since TMZ metabolism is highly pH dependent, we moni- 
tored the medium pH (pH 0Ut ) and found an initial increase 
followed by a return close to its initial value which was mod- 
eled as (Figure 3g): 

P H out (f) = pH 0+ e-^-e- c ^ 

in which pH 0 is the initial pH value, and £ pH and C pH are two 
phenomenological parameters. Available experimental data 
in U87-cultured cells from 20 allowed us to derive the following 
formula for the intracellular pH (pH in ) (see Supplementary 
Data): pH in = pH 0Ut -0.01 . 



Table 1 Model parameters 



Parameters (unit) 



Reaction 



Buffer 
solutions 



Cell 
culture 



Mice 



Humans 



Metabolism 

fcr„(h-') 

\ (h- 1 ) 

TMZ transport 

p T , pj umor (l.h- 1 ) 
p T2 , p™ (l.h- 1 ) 

pbrain (| h ~1) 

p T f in (l.h" 1 ) 

gtumor (|.h- 1 ) 

g} u 2 mor (l.h- 1 ) 

g, brain (| 

ctff 1 (l.h- 1 ) 
in vitro experiments 
TMZ Q (|jmol/l) 
P A (l-h- 1 ) 
P A2 (l-h- 1 ) 
PH 0 

B pH (h- 1 ) 

C p h (h") 
In vivo blood forcing 
C 0 (Mmol/l) 
K a (h- 1 ) 

*c,ea r ('-h- 1 ) 

Pharmacodynamics 
Fixed parameters 

*ca,(h- 1 ) 



TMZ metabolism into MTIC 
MTIC degradation into AlC 



TMZ brain tumor influx 

TMZ brain tumor efflux 

TMZ normal brain influx 

TMZ normal brain efflux 

TMZ blood to tumor interstitial fluid transport 

TMZ tumor interstitial fluid to blood transport 

TMZ blood to normal brain interstitial fluid transport 

TMZ normal brain interstitial fluid to blood transport 

Initial TMZ concentration 
AlC U87 cell influx 
AlC U87 cell efflux 
Initial medium pH 
pH parameter 
pH parameter 
function 

Initial TMZ gut concentration 
First-order absorption 
First-order elimination 

DNA adducts formation rate 



1.1x10-' 
1.96 
292 
0.33 



2.09 (0.02) 
0.31 (0.26) 



0.0038 (16.9) 0.039(16.9) 
0.008 (17) 0.081 (17) 
0.127(16.9) 
0.154 (51.9) 
0.21 x10- 3 (22.4) 
0.15x10- 3 (47.4) 
0.24 x10- 3 (29.2) 
0.49x10- 3 (47.4) 

59.2 (0.04) 
0.0034 (28.5) 
0.0072 (28.5) 
7.57 (0.02) 
14.6 (0.56) 
998.3 (1 .4) 

2330.6 (5.5) 
0.64 (28.1) 
18.7(30.5) 



1.81 (0.33) 



Naive Estimated 
17.5 (16.9) 
36.32 (17) 

507.2 (16.9) 

620.3 (51.9) 

0.094(22.4) 0.019 (122) 

0.067(47.4) 0.069(170) 

0.97(29.2) 0.197(122) 

1.97(47.4) 2.004(170) 



Methylating cation degradation 



6,000 



334.6(107.2) 
0.398 (35.2) 
2.73 (99.2) 



Mean values (%CV) are shown. Bold numbers indicate that parameters were scaled from the 
AlC, 4-amino-5-imidazole-carboxamide; MTIC, metabolite 5-(3-methyltriazen-1-yl)imidazole-4 



previous scale. 

-carboxamide; TMZ, temozolomide. 
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Figure 2 TMZ pH-dependent conversion in buffer solutions, (a-d) TMZ, MTIC, and AIC concentration-time profiles in buffer solutions at 
indicated pH. Solid lines and dots correspond respectively to the best-fit profiles and data from ref. 15 of TMZ (black), MTIC (orange), and AIC 
(yellow), (e) Total methylating cation, defined as the sum of the cation C and all products of its degradation, with respect to buffer solution pH 
for different exposure duration. AIC, 4-amino-5-imidazole-carboxamide; MTIC, metabolite 5-(3-methyltriazen-1-yl)imidazole-4-carboxamide; 
TMZ, temozolomide. 



The final model (Figure 1b, equations in Supplemen- 
tary Data) required 14 parameters to be determined: 4 
parameters for TMZ pH-dependent metabolism (k JQ , X v 
*W K)> 4 transport rate constants (p T p T2 , p A p A2 ), 3 pH 
parameters (pH 0 , £ pH , C pH ), the degradation rate of the 
cation /c cat ,1 PD parameter (/c add ) and the initial TMZ extra- 
cellular concentration from the U87 cell studies (TMZ Q ). 
Parameters k caV k J0 and k m were directly inferred from 
the buffer study whereas X T and A, M were estimated for the 
in vitro system allowing a 50% deviation from the buffer 
solution values. The 1 1 remaining parameters were esti- 
mated by fitting TMZ, MTIC, and AIC extra- and intracellu- 
lar concentration-time data and extracellular pH results in 
U87 cells together with available PD data from which it was 
estimated that an exposure of TMZ at 800 umol/l during 



1 h led to a DNA adduct concentration of 70 nmol/l in the 
absence of DNA repair mechanisms 21 (see Supplementary 
Data). 

Agreement between fitted profiles and observed values 
together with low %CV of parameter estimates support model 
validity (Figure 3; Table 1). TMZ exposure led to increasing 
concentrations of AIC. DNA adduct concentration was pre- 
dicted to increase as no repair mechanisms were considered 
in this model. 

TMZ brain disposition in mice 

The next step consisted of designing a mechanistic PK 
model of TMZ in brain, based on two complementary inves- 
tigations; a microdialysis study in normal mice, which pro- 
vided serial measurements of TMZ plasma and normal brain 
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Figure 3 TMZ PK in U87 cell culture, (a-f) TMZ, MTIC, and AIC extra- and intracellular concentration-time profiles. Dots represent 
experimental data (mean ± SD); solid lines are the mathematical model best fit. (g) Medium pH time evolution during TMZ exposure. 
Experimental data points are the mean and SD of three independent experiments, (h) DNA adducts concentration simulated by the trained 
model. AIC, 4-amino-5-imidazole-carboxamide; MTIC, metabolite 5-(3-methyltriazen-1-yl)imidazole-4-carboxamide;TMZ, temozolomide. 



interstitial fluid concentrations, and a steady-state investi- 
gation in which TMZ plasma, whole brain, and brain tumor 
concentrations were measured in mice bearing intracerebral 
U87 tumors. 

Because we focused on TMZ brain disposition, its blood 
concentration T b , 00d was accounted for by a forcing func- 
tion (Figure 1c) chosen as the best-fit to the plasma mea- 
surements from the microdialysis mouse study assuming a 
blood:plasmaTMZ concentration ratio of 1 2223 : 



d7" b | ( 



--C 0 K a e~ 



^clear ' blood 



where C 0 , K , and k . are the initial TMZ gut concentration, 

0' a' clear ^ ' 

the first-order absorption rate constant, and the first-order 
elimination rate constant, respectively. C 0 was computed 
for each individual mouse as the total administered dose 



divided by the gut volume, assumed equal to 1.7 ml, 24 
whereas K a and /c clear were estimated from plasma data 
(Table 1; Figure 4a). 

The normal brain and brain tumor were assumed to be con- 
stituted of homogeneous normal and U87 glioma cell popula- 
tions, respectively, and were represented by physiologically 
based PK models, which consisted of interstitial and intra- 
cellular compartments of volume \4 b [ e a r n , ^ n b [ r a a in , K£ or , and 
V^ntra or (Figure 4a). Utilizing the model designed in the buf- 
fer solution study, TMZ pH-dependent metabolism into MTIC 
was represented in all brain compartments, whereas MTIC 
and the methylating agent — which were assumed unable to 
diffuse from one compartment to another — were only explic- 
itly modeled in the intracellular compartments in which the 
methylating cation creates DNA adducts at the rate /e add . The 

A,., and the PD rate 

M 



metabolic parameters /c cat , k J0 , X v k M 
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Figure 4 TMZ PK in mice, (a) TMZ blood concentration-time profile. Solid line is the forcing function best fit. Open circles are the measured 
TMZ plasma concentrations from the microdialysis study, (b-e) TMZ interstitial (b) and intracellular (c), intracellular methylating cation 
(d), and DNA adducts (e) concentration-time profiles in normal brain (solid line) and brain tumor (dashed lines). Open circles (b) are TMZ 
normal brain concentrations from the microdialysis study, (f) TMZ concentrations in the whole normal brain and brain tumor after steady- 
state dosing regimen. Open circles are the experimental results in brain tumor-bearing mice; black diamonds are the model best fit. PK, 
pharmacokinetics; TMZ, temozolomide. 



/c add were directly inferred from the U87 cell study. pH values 
were assumed identical in the normal brain and tumor region, 
as experimentally demonstrated in humans 25 , and were set 
to 7.2 in the intracellular compartment, 7.3 in the interstitial 
fluid, and 7.4 in the blood. 2426 

The tumor cell membrane transport parameters p| umor and 
PT 2 mor were computed by scaling the in vitro parameters p T 
and p T2 using the following formulas, which suppose that 
the contact surface between the cells and the extracellular 
medium was proportional to the total cell volume: 



■-.tumor _ „ 

Pt - Pt ' 



V 'tumor 
intra 



and p| u 2 mor = p T2 



i/tumor 
^intra 



in which V ]n is the in vitro intracellular volume, and V; n tra or 
the tumor intracellular volume computed as a fraction of the 
total tumor volume: Vj»™ r = Ua ■ KotaT ■ 

The six remaining parameters qT rain , qjf n , pj a ' n , Pji' n , 
q { j mor , andqT 2 mor were estimated by utilizing the micro- 
dialysis and steady-state mouse studies. Since these 
datasets were not sufficient to ensure unique parameter 
estimates, two constraints were added assuming that drug 
transport from interstitial fluid to intracellular medium, and 
from interstitial fluid to blood, were equivalent in the brain 

tumor and normal brain such that: n brain _ j4rajn m n tumor and 
V Pt - Pt 

^brain _ "brain stumor "tumor 
WT2 - — 4T2 ' 



While TMZ microdialysis normal brain interstitial fluid 
measurements were directly compared to the correspond- 
ing model variable, whole normal brain and brain tumor 
TMZ concentrations did not explicitly appear in the model 
and were computed as follows. We assumed that the normal 
brain and tumor were exclusively constituted of blood, inter- 
stitial fluid, and intracellular components with respective vol- 
ume fractions /biood, /i nter , and /j ntra set to 2, 20, and 78% for 
the normal brain and 5, 30, and 65% for the tumor. 27 " 29 TMZ 
concentrations in the whole normal brain or tumor can then 
be computed as (omitting superscripts for normal brain and 
tumor): r tota , = /bioodTbiood + Winter + Tintra ■ During steady- 
state administration, TMZ concentrations in brain compart- 
ments rapidly reached quasi-steady state, as suggested 
by the data and adjoining simulations, and were computed 
S-TSood with A = -^— and 

PT2 + K T 



as: T* tra 
B = 



AT* T* 

' ^ ' inter > ' inter 



q T2 +/c| nter +p T 



^--PT2 

V. V. 
v inter v inter 



Kntra A where T* ood Was Set 



to the steady-state plasma measurement in each mouse. 
Volumes Vj nter and V Mra were computed for the normal brain 
and the tumor as V[ nter = ^ter^otai and \4 tra = Ua^totah the 
whole normal brain volume being set to 360 ul, 24 and the 
tumor volume to 111.5 ul, its mean value in the steady- 
state study. Then, TMZ steady-state total concentration 
was computed as: 7^ = (/ b , 00d + f Mer A + U a AB) ■ T^ 00d ■ 
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Figure 5 TMZ PK in brain tumor patients, (a) TMZ blood concentration-time profile. Solid line is the forcing function best fit. Open circles areTMZ 
plasma concentration measurements in patients. 31 (b-e) TMZ interstitial (b) and intracellular (c) concentration-time profiles in the normal brain 
(solid lines) and brain tumor (dashed lines). Blue color corresponds to the model calibrated with naive parameters, whereas black color represents 
the best-fit model with estimated parameters. In (b), open diamonds areTMZ normal brain concentrations from microdialysis measurements 
in patients, 31 which are shown in greater details in the inset, together with TMZ normal brain interstitial concentration best-fit profile, (f) Global 
analysis of parameter sensitivity on drug efficacy in brain tumor (see Results and Methods). PK, pharmacokinetics; TMZ, temozolomide. 



The parameter estimation procedure consisted of two 
independent steps leading to the estimation of, first, qT rain , 
q^ ain , and p^ ain and secondly gf™ (Table 1). Normal brain 
parameters were estimated by combining TMZ whole nor- 
mal brain and brain microdialysis measurements as follows. 
Since the steady-state and microdialysis studies were per- 
formed in 7 and 10 mice respectively, a parameter set was 
estimated for each of the 70 possible pairs composed of one 
individual of each study, by fitting the individual steady-state 
normal brain and plasma TMZ concentrations together with 
the individual normal brain interstitial TMZ concentrations, 
setting the blood forcing function parameters to their values 
estimated for the corresponding mouse of the microdialy- 
sis investigation. The resultant 70 parameter sets were then 
used to calculate the mean and %CV of each parameter. 
Finally, q| umor was independently estimated by fitting indi- 
vidual steady-state plasma and whole tumor concentrations, 
and computing the mean and %CV of individual parameters. 
The model best fit and experimental data agreed (Figure 4). 

Both in the normal brain and tumor, transport parameters 
between interstitial and intracellular compartments were 
three orders of magnitude greater than the ones between 
blood and interstitial fluid (Table 1). Intracellular transport 
was faster in the normal brain than in the tumor, however 

n tumor n brain 

parameter ratios ^T2 and ^T2 were in the same range. 

stumor ~>brain 

Pi Pt 



stumor 
4T2 
~ tumor 



was 2.5-fold lower than 



^brain 
HT2 
^brain 
HI 



which was 



consistent with BBB disruption in the tumor often reported 
in the literature. 30 As a consequence, best-fit TMZ intersti- 
tial and intracellular concentrations were nearly equivalent in 
the normal brain with an area under the curve (AUC) ratio of 
1 .2, whereas, in the tumor, TMZ interstitial concentration was 
greater than the intracellular concentration as reflected by an 
AUC ratio equal to 2.08. Finally, the TMZ intracellular AUC 
was 1.7-fold greater in the brain tumor than in the nor- 
mal brain which led to a 1.8-fold higher final DNA adduct 
concentration (Figure 4). 

TMZ PK in brain tumor patients 

The final step of our multiscale approach was directed at 
predicting TMZ brain disposition in brain tumor patients. For 
the human study, the structure of the mouse model was kept 
identical and parameters were adapted as follows. First, the 
blood forcing function parameters were determined from 
mean and standard deviation values of TMZ plasma concen- 
trations available in seven brain tumor patients after a single 
oral dose of 150mg/m 2 . 31 Parameters of TMZ metabolism 

(^caf ^TO' K *W and K) and PD (*add)' 3S We " 3S P H Va ' UeS 

and brain and tumor volume fractions were inferred from the 
mouse study. Volumes of whole normal brain and brain tumor 
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Figure 6 A multiscale approach to design cell-type-specific PK-PD 
models for personalized medicine. The in vitro study is based on 
experimental drug investigations in solutions and in cell culture, 
which allow the design and calibration of a mathematical model 
describing the extracellular and intracellular drug PK-PD. Next, for 
the preclinical scale, the model of the interstitial and intracellular 
compartments is directly inferred from the in vitro study and 
parameters are either kept identical, scaled, or estimated utilizing 
whole animal or tissue-based investigations that in our case included 
two separate studies: normal brain microdialysis and brain tumor 
orthotopic studies. This step allows the design and calibration of 
a whole-body or tissue model that designate specific tissues and 
cells of interest. For humans, the whole-body model or specific 
tissue structure remains and parameters are either kept identical, 



were set to 1 .45 I 24 and 0.05 I, 32 respectively. Cell membrane 
transport parameters were scaled up from mice to humans 
by assuming that the contact surface between interstitial and 
intracellular mediums was proportional to the cell volume 
such that: 

i/human 
phuman _ pimouse # "intra 



V 'mouse 
intra 

where p stands for p£ rain , p^ ain , p| umor , or p| u 2 mor , and V[ ntra 
for the normal brain or tumor intracellular volume. 

Transport parameters from blood to interstitial fluid were 
either directly scaled up from the mouse study — referred to 
as the naive estimation — or estimated using TMZ normal 
brain concentrations measured by intracerebral microdialy- 
sis in the same seven patients from whom plasma measure- 
ments were available. 31 For the naive parameter scale-up, the 
following formula was used, which assumes that the contact 
surface between the blood and the interstitial fluid was pro- 
portional to the volume of the blood vessels: 



--q" 



i/human 
"blood 



g£? in , or q| u 2 mor and V Hl 



where q stands for gf £ , , ui anu v b | 00d 

for the volume of the blood compartment in the nor- 
mal brain or tumor. For the parameter estimation based 
on human data, qj ain and qj 2 ain were computed by fit- 
ting the model to the brain microdialysis TMZ concentra- 
tions from patients. Tumor parameters were subsequently 
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deduced by assuming, that q*f 

mouse study and that the ratio between brain and tumor 
BBB permeability was identical in mice and in humans: 

stumor i/brain i/tumor 
stumor _ n brain WT,mouse "total,mouse "total, human /Tahlp "H 
WT.human WT.human brain ./tumor i/brain \icxuiKSi). 

4T,mouse "total, mouse "total, human 

Naive scale-up overestimated TMZ concentrations in the 
normal brain and tumor as AUC values of best-fit TMZ inter- 
stitial and intracellular concentrations were on average 5.1- 
fold greater than the ones of the human data-fitted model 
(Figure 5a-e). Concerning the human data-fitted model, the 
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ratio tumor was 2.8-fold lower than ~bra¥T which suggested a 

disruption of the BBB in the tumor region as in the mouse 
study. Interestingly, the ratios between TMZ interstitial and 
intracellular AUC were the same in the preclinical and human 
data-fitted simulations. As a consequence of having inferred 
from the mouse study the ratios between brain and tumor 
transport parameters, TMZ intracellular AUC ratio between 
brain and tumor was close to the preclinical one as it was 
equal to 1.7, leading to a 1.7-fold higher final DNA adduct 
concentration in the tumor compared to the normal brain 
(Figure 5). The %CV of parameters were relatively high due 
to the inherent variability of the clinical data that was uti- 
lized for parameter estimation (Figure 5a, b; Table 1). Next, 
we performed a global sensitivity analysis to determine the 

scaled from the preclinical study, or estimated from clinical data 
in order to compute an average population parameter set. Finally, 
the human model is partially re-parameterized from available data 
on the individual patient and utilized in optimization procedures, 
which provide personalized drug combinations and administration 
schemes. PD, pharmacodynamics; PK, pharmacokinetics. 
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important parameter regarding the efficacy of TMZ computed 
as the AUC of DNA adduct concentration in the tumor. The 
parameter g^ 01 " which corresponds to the BBB had the 
highest influence followed by /c clear , TMZ clearance rate in the 
plasma (Figure 5f). 

DISCUSSION 

In this study, we undertook a multiscale approach to quantita- 
tively characterize TMZ brain disposition in patients. Physio- 
logically based mechanistic models were used to describe at 
the molecular scale the biochemical events determining TMZ 
PK. The main challenge of this type of modeling lays in the 
estimation of the model parameters which may not be directly 
measurable from clinical investigations. However, since those 
parameters do have a physical meaning and correspond to 
biochemical reactions, one may evaluate them independently 
in in vitro and in vivo studies and subsequently scale them to 
patients (Figure 6). Once population parameters are deter- 
mined for humans, the model can be partially reparameter- 
ized for a specific patient allowing the design of personalized 
drug combinations and administration scheme. 

To arrive at the final TMZ brain disposition model, a pipeline 
of in vitro and in vivo investigations were needed that con- 
sisted of four main scales; in vitro aqueous solution and U87 
glioma cell studies to characterize membrane transport and 
metabolism, in vivo mouse studies that allowed designing and 
calibrating the brain disposition model, and finally the clinical 
investigation for which the mouse model structure was main- 
tained yet revised by parameter estimation. At each step, the fit- 
ted models closely agreed with the observed data and yielded 
reliable parameter estimates (low %CV) which supports the 
use of this multiscale strategy. It should be appreciated that 
the sole use of either the microdialysis or steady-state mouse 
studies would not have provided reliable parameter values, 
and without the in vitro studies, TMZ membrane transport and 
metabolism would not have been accurately estimated. 

The current TMZ brain disposition model, whose tumor 
intracellular compartment represents a homogenous popula- 
tion of glioma cells, can be extended by considering several 
intracellular compartments within the tumor, thus represent- 
ing different cell types, that may have important roles in tumor 
progression, therapeutic targets, and drug resistance. 8 Fur- 
ther, models of organs of toxicity, again with subpopulations 
of cells of different susceptibilities, could be used in conjunc- 
tion with the cell-type-specific tumor models to refine drug 
treatment regimens. This can readily be appreciated for TMZ 
whose main dose-limiting toxicity of myelosuppression could 
be considered by adding to the model bone marrow compart- 
ments representing hematopoietic stem and progenitor cells. 
A hypothetical representation of such model is provided in 
Figure 1d in which the blood and interstitial fluid compart- 
ments drive the dynamics of different intracellular compart- 
ments; highlighted in tumor and bone marrow. This type of 
intracellular PK-PD models including both target and toxicity 
organs provides a quantitative tool to evaluate therapeutic 
indices in a mechanistic manner. 

It is appreciated that deriving a multiple cell-type model 
will require in vitro investigations — as done here for glioma 
cells — to delineate PK-PD behavior in cell types of interest. 



Apart from the direct use of in vitro cell studies, isolation of 
different cell fractions from ex vivo samples may be afforded 
by cell sorting. 33 While drug concentrations may be altered 
by cell sorting procedures, certain PK-PD features especially 
slowly reversible or irreversible endpoints such as apopto- 
sis, may not be unduly altered. Sorting untreated xenograft 
tumors into specific cell fractions provides an alternate means 
to initiate in vitro PK-PD studies, and when coupled to brain 
microdialysis and whole tissue studies as done here provide 
the means to obtain cell-type-specific models. Approaches 
to dissect bone marrow into hematopoietic stem and progeni- 
tor cells could include collection of bone marrow from mice or 
collection of ex vivo patient samples that can be cultured for 
further drug studies. 3435 Thus, the development of cell-type- 
specific PK-PD models is feasible and likely to benefit from 
further technical innovations to enable multiscale strategies. 

Personalized medicine and inherent tumor heterogene- 
ity provide strong motivation to pursue the development of 
cell-type-specific intracellular PK-PD models. The proposed 
pipeline of investigation is offered here as an essential first 
step to achieve this goal that should contribute to under- 
standing drug action among heterogeneous cell types. 

METHODS 

In vitro cellular PK studies of TMZ 

U87 human glioma cells (ATCC, Manassas, VA) were cultured 
in Dulbecco's modified Eagle's medium supplemented with 
10% fetal bovine serum, 100 units/ml penicillin, and 100 units/ 
ml streptomycin and maintained in a humidified atmosphere of 
5% C0 2 in air at 37 °C. For PK experiments, 10 6 cells were 
seeded onto 60 mm Petri dishes, cultured overnight and TMZ 
(50 umol/l) added, and the dishes sampled for 2 h. At each time, 
dishes were placed on ice instantly, and 1 ml of medium was 
collected and kept on ice while the reminder of the medium was 
aspirated. The cells were washed twice with 1 ml of cold phos- 
phate-buffered saline, scraped in 200 ul of cold MeOH and kept 
on ice. For the detection of TMZ which is stable at pH < 4, 1 0 ul 
of medium was transferred into 300 ul of acid internal standard 
solution (0.1 ug/ml dacarbazine and 0.1% formic acid in MeOH) 
and 90 ul of MeOH. Another 10 ul of medium was transferred 
into 300 ul of internal standard solution without formic acid and 
90 ul of MeOH for MTIC and AIC detection which are more 
stable at alkaline pH.The cells were sonicated (6 s, three times) 
and 20 ul of cell lysate were added to 60 ul of internal stan- 
dard solution with or without formic acid as for medium samples 
and kept on dry ice. For pH measurements, cells were incu- 
bated with TMZ (50 umol/l) and the medium pH was instantly 
measured. All medium and cell lysate samples for TMZ, MTIC, 
and AIC were analyzed on the same day by high-performance 
liquid chromatography-tandem mass spectrometry (LC/MS/ 
MS) using the multiple reaction monitoring mode specific to 
each analyte and the internal standard. The LC/MS/MS sys- 
tem (described in Supplementary Tables S1-S3 online) was 
optimized for each compound in positive ionization mode. The 
analytical methods were specific and sensitive with a lower limit 
of quantification of 0.01 umol/l. The intraday and interday vari- 
abilities were less than 15% in all matrixes. The average run 
time was ~5min. Raw mass spectrometry measurements were 
normalized to 1 million of cells (see Supplementary Data). 
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Microdialysis experiments in nude mice 

Four days before mice entered the PK study, a normal brain 
microdialysis guide cannula was implanted similar to that previ- 
ously described for rats. 32 One day before the PK sampling, a 
catheter was inserted into the right carotid artery of each mouse 
under isoflurane anesthesia for blood sampling. On the day of 
and prior to TMZ administration, the dummy cannulas were 
replaced by microdialysis probes (MRB-1-5; Bioanalytical Sys- 
tems, West Lafayetter, IN) and the probe TMZ recovery values 
were determined using the retrodialysis method as described 
previously with minor modifications. 32 The flow rate was main- 
tained at 1 ul/min during the washout and microdialysate frac- 
tions were automatically collected into individual 250-ul plastic 
vials containing 5 ul of 0.1% formic acid every 20min with a 
CMA/170 refrigerated fraction collector. The percent relative 
recovery of TMZ from the normal brain was 8.1 ± 1 .0%. Animals 
were given a single oral dose of 20mg/kg of TMZ with blood 
samples (15 ul) collected over 6h. Plasma was separated by 
centrifugation and stored at -80 °C until analysis. TMZ concen- 
trations in plasma and normal brain microdialysate samples were 
determined using a validated reversed-phase high-performance 
liquid chromatography method with ultraviolet detection. 36 

Steady-state experiments in nude mice 

A previously described orthotopic mouse model of human 
glioma was used with minor modifications. 37 Briefly, anesthe- 
tized NIH-Swiss nude mice (nulnu) were secured in a stereo- 
taxic instrument, and a 1.5cm longitudinal incision was made 
to expose the skull. Using a sterile 10-ul Hamilton syringe with 
a 26-gauge needle attached to the stereotaxic frame, U87MG 
cells (1 x10 6 ) in 10 ul of saline were injected 2.5 mm deep into 
the left caudate putamen over a 5-min period. The skin was 
sutured to close the wound. After surgery, animals were returned 
to their cages and received a regular diet and water ad libitum. 

Twenty-four days after tumor cell implantation, tumor- 
bearing mice were prepared for aseptic surgery to introduce 
a right carotid artery cannula for intraarterial infusion of TMZ 
(dissolved in 0.9% NaCI containing 25% dimethyl sulfoxide) 
administered the next day over 1 80 min at a rate of 1 33.3 ug/kg/ 
min preceded by a tail veil injection of 7.6mg/kg of TMZ; both 
designed to achieve steady-state plasma concentration of -17 
ug/ml. 37 At the end of the infusion, mice were anesthetized and 
blood samples were taken from the vena cava. Plasma was 
prepared by centrifugation of the heparinized blood and stored 
at -80 °C. After the animals were sacrificed, normal brain tis- 
sues and brain tumors were immediately excised, snap frozen 
on dry ice and stored at -80 °C.TMZ concentrations in plasma, 
normal brain, and brain tumor homogenates were determined 
using a validated reversed-phase high-performance liquid 
chromatography method with ultraviolet detection. 36 

Model development, parameter estimation, and sensitiv- 
ity analysis 

TMZ PK models are based on ordinary differential equations, 
which were solved by the Matlab function odes15 (Math- 
works, Natick, MA). Chemical reactions were modeled by 
the law of mass action and passive drug transport by Fick's 
first law of diffusion. 3839 Parameter estimation was performed 
using a weighted least square approach in which the cost 
function was minimized by the Covariance Matrix Adaptation 



Evolution Strategy algorithm. 4041 For the U87 and patient 
studies, Monte Carlo simulations were utilized to compute 
parameter estimates and %CV (see Supplementary Data). 

For the global sensitivity analysis of the human model, all 
parameters' lower and upper bounds were set to 100-fold 
lesser and greater than their estimated values and 18,000 
parameter sets were generated from cross-sampling by Salt- 
er's extension of Sobol's method using the MOEA frame- 
work (version 2.0). 42 For each parameter set, the AUC of 
DNA adduct concentration in the tumor during the 1 8 h follow- 
ing an oral dose of TMZ of 150mg/m 2 was calculated using 
Matlab. Finally, parameters' total-order sensitivity indices and 
their confidence intervals were computed using the Sobol 
analysis of the MOEA framework. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

Although it has been experimentally demonstrat- 
ed that heterogeneous cell populations both within 
the brain tumor and the organs of toxicity display 
different drug responses, there has not been any 
modeling effort to optimize temozolomide (TMZ) 
chemotherapy on a cell-type-specific basis. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

y This study aims at designing a physiologically 
based model of TMZ brain disposition through 
a multiscale modeling approach, thus provid- 
ing a critical tool toward the personalization of 
TMZ-based treatments. 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y This study presents the first mechanistic model 
of TMZ brain disposition in mice and in brain 
cancer patients, and provides a reliable basis to 
design cell-type-specific models of TMZ PK-PD. 



HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y Cell-type-specific models of TMZ PK-PD may 
be used as comprehensive tools to personalize 
TMZ-based chemotherapies allowing the de- 
sign of optimal drug combinations and adminis- 
tration schedules for individual patients. 
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